Does strange kinetics imply unusual thermodynamics? 
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We introduce a fractional Fokker-Planck equation (FFPE) 
for Levy nights in the presence of an external field. The equa- 
tion is derived within the framework of the subordination of 
random processes which leads to Levy flights. It is shown 
that the coexistence of anomalous transport and a potential 
displays a regular exponential relaxation towards the Boltz- 
mann equilibrium distribution. The properties of the Levy- 
flight FFPE derived here are compared with earlier findings 
for subdiffusive FFPE. The latter is characterized by a non- 
exponential Mittag-Leffler relaxation to the Boltzmann dis- 
tribution. In both cases, which describe strange kinetics, the 
Boltzmann equilibrium is reached and modifications of the 
Boltzmann thermodynamics are not required. 
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Strange kinetics which involves diffusional anoma- 
lies, both sub- and superlinear, and nonexponential re- 
laxations is quite wide-spread and has been observed in a 
broad range of systems [0J^~|^] ■ The ubiquity of strange 
kinetics rests upon generalization of the central limit the- 
orem due to Levy 0, a generalization that puts heavy- 
tailed distributions on a same level of importance as the 
well-known Gaussian. 

Anomalous diffusion in the presence or absence of an 
external field has been modelled in a number of ways, 
including fractional Brownian motion generalized 
diffusion equations g], continuous time random walk 
(CTRW) models jl(J , Langevin and generalized Langevin 
equations jnj and generalized thermostatics jl2j . In par- 
ticular CTRW has been demonstrated to be a powerful 
approach in describing subdiffusive as well as superdif- 
fusive processes and in interpreting experimental results, 
ft is not straightforward, however, to incorporate force 
fields and boundary conditions in this formalism. 

An alternative approach to processes which display 
strange kinetics is based on fractional equations, which 
are suitable for handling external fields and for consider- 
ing boundary value problems. In the case of subdiffusion 
it was realized that the replacement of the local time 
derivative in the diffusion equation by a fractional oper- 
ator accounts for memory effects responsible for anoma- 
lous behavior |]|n|. In the presence of an external field 
a fractional Fokker-Planck equation FFPE has been in- 
troduced [EUil: 
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o-D t " is a fractional Riemann-Liouville operator < 
a < 1, and Km is a generalized (sub-)diffusion coeffi- 
cient, having the dimension [-KV a )l = [L 2 /t Q ] . The force 
f(x) is related to the external potential U(x) through 
f(x) = —dU/dx, and ks is the Boltzmann constant. The 
differential operator aD\~ a acting on functions of time 
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The FFPE, Eq.(0), has been derived using a Kramers- 
Moyal expansion of the CTRW nonlocal equation [jUg . 
The solution of this FFPE is characterized by a subdif- 
fusive behavior and by a nonexponential Mittag-Leffler 
decay of the single modes. The decoupled structure of 
Eq. ([!]) guarantees that the Boltzmann distribution is at- 
tained at equilibrium H|l3|T]|]. We note that the latter 
is also a property of the regular Fokker-Planck equation 
corresponding to a = 1. 

Less clear has been the situation for FFPEs which 
correspond to Levy spatial flights. Previously proposed 
equations |2"HIl{| seem not to lead to the Boltzmann distri- 
bution, a point whose impact has been overlooked. This 
might suggest therefore that strange kinetics requires un- 
usual thermodynamics |l2|. Here we derive a FFPE for 
Levy flights in the presence of an external force. Our 
starting point is a representation of Levy flights in terms 
of a subordination of random processes |L5 16 1. This 



where Cfp is the Fokker-Planck operator, 



representation corresponds to processes in which space 
and time are decoupled and it does not account for Levy 
walks [p]PJTo[|. Namely, in what follows we obtain a di- 
verging mean-square displacement in the force-free case. 
The solution of the FFPE which we derive again leads to 
the Boltzmann distribution in the equilibrium limit, re- 
emphasizing that there is no need to modify conventional 
thermodynamics in order to obtain strange kinetics. We 
bring some examples for solving this FFPE for boundary 
value problems. 

As we proceed to show, the corresponding generaliza- 
tion of the Fokker-Planck equation for Levy flights is: 

^P(x, t) = -JfW {-C FP ) a P(x, t), (4) 

where the operator (—Cfp) is the a-th power of the 
operator —Cfp = —d 2 /dx 2 + d / dx(f [x) / UbT) as will 
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be derived below, and the corresponding generalized 
(super-)diffusion coefficient has as dimension [iff")] = 
[L 2 «/t]- 

The CTRWs can be viewed as Markovian random 
walks on a lattice (with lattice constant a) given in terms 
of the number of steps n of the random walker. P(x, n) 
is a probability distribution function (pdf) of the parti- 
cles' displacement x after n steps. The number of steps n 
performed during the time t follows the probability dis- 
tribution S(n, t), which may include memory effects [ [l7|| . 
The overall displacement during time t is then given by 



P(x,t) 



^2P(x,n)S(n,t). 
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In the force-free case the pdf P(x,n) corresponds typi- 
cally to normal diffusion behavior, and thus x 2 oc n. On 
the other hand, the typical number of steps can grow sub- 
or superlinearly in time, so that the overall behavior can 
be anomalous. 

Here we concentrate on the superdiffusive case and as- 
sume that the random process {n(t)} is characterized 
by a diverging mean density of events, so that the first 
moment of the number n of steps does not exist. As a re- 
alization of such a process we can take that the numbers 
of jumps during different time intervals of unit length 
are independent random variables distributed according 
to S(n, 1) oc n _1_a . For t large enough the distribution 
S(n,t) tends then to a stable Levy-law L(n;a,/3) [ p"5| . 
Since n is nonnegative, this law is the one-sided extreme 
distribution for which (3 = —a (0 < a < 1). If different 
time intervals t are considered, the distribution S(n, t) 
scales as 
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Imagine now a random walker moving under the influ- 
ence of a weak force f(x). Such a force introduces an 
asymmetry into the walker's motion, since the probabil- 
ities for forward and backward jumps, w + and W- are 
now weighed with the corresponding Boltzmann-factors, 
w + /w- = exp(fa/k B T). For small / one can take 
w + = 1/2 + fa/2k B T and w_ = 1/2 - fa/2k B T. Note 
that the process described in such a way is a Markovian 
one, and can be characterized by a transition probability 
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P(x-x',n)S(n,At). (7) 



For At in the intermediate range, i.e. large enough 
to view both x and n as being continuous and 
to approximate P(x,n) by the Gaussian P(x,n) — 
(27m)~ 1/2 exp(-(x-un) 2 /2a 2 n) with v = fa/2k B T, 
yet small enough to have the typical displacement small 
on the scale of change of f(x), one obtains: 

W(x,t + At\x',t) = 
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S(n, At)dn. (8) 



The overall Markovian process is then governed by the 
integral Chapman-Kolmogorov equation 

Let us concentrate on the long-time, large x, behavior of 
the system and take the force / to be smooth. On such 
scales we consider the x- Fourier-transform of Eq. (^) and 
get: 

P(k,t + At) = (10) 

i 

dnexp (-(ikv + k 2 )a 2 n) S(n, At)P(k', t). (11) 

For random processes leading to diffusive behavior, the 
first moment of the distribution S(n, At) for small At 
exists, so that one can expand for small k the exponential 
into a power series, getting: 

/>oo 

P(k,t + At)= / dn(l- (ikv + k 2 )a 2 n) S(n,At)P(k,t) 
Jo 
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= P(k, t) - (ikv + k 2 )a 2 (n) At P(k, t). (12) 

(a Kramers- Moyal procedure). For normal diffusive pro- 
cesses one has (n) At ~ wAt, where w is the jumping rate, 
so that in the continuum limit 

^P(k,t) = -K^k^ + k 2 \p(k,t), (13) 

with K = a 2 w/2 being the diffusion coefficient. In the 
x-representation this is the conventional Fokker-Planck 
equation (FPE) |J): 
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In the case when S(n, At) is a Levy-stable law of in- 
dex a, < a < 1, the first moment of n diverges, 
and the series expansion of the exponential, Eq.(|l2|), is 
not possible. On the other hand for a < 1 the inte- 
gral 4>(k) — L exp(— kt)S(t, At)dr converges for each 
K = £ + if)i Re£ > 0, and is a stretched-exponential 
function . For extreme Levy-stable distributions with 
< a < 1 (those which vanishing identically for negative 
arguments) one has (j)(k) — exp(— n a ). Thus, performing 
the integration in Eq.(|10|) one gets: 



P(k, t + At) = exp 
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Expanding now the exponential and repeating the steps 
leading to Eq.^) we have 
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+ k 2 ) P(k,t) (16) 



Comparing the terms —(ikf/ksT + k 2 ) a and 
— (ikf/ksT + k 2 ) in Eqs. jlq ) and (|l^), which represent 
the corresponding transport operators £ Q and Lfp = £1 
in Fourier space, we see that they are connected by the 
relation C a = — (— Cfp)" ■ The same relation holds, of 
course, when one shifts to the ^-representation: 
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see note Jl9| ] . Note that Eq.(|l7j) differs from the expres- 
sions proposed in p"i"| ] where either only the second part of 
the Fokker-Planck operator (a A-term) is changed (and 
corresponds in our notation to — (—d 2 /dx 2 ) a ) or where a 
sum of two terms is introduced, so that fractional space 
derivatives of the orders a and 2a appear. Note that 
in general C a can not be decoupled into additive parts 
responsible separately for drift and for diffusion. 

Some important properties of the Levy-diffusion in the 
presence of a force field stem from Eq. (p"7j) . Since —Lfp 
and C a commute with each of their powers, the eigen- 
functions of these operators coincide. The corresponding 
eigenvalues of C a are those of —Cfp raised to the power 
of a: 
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Note that the eigenfunctions of — (ikf/ksT + k 2 ) and of 
— (ikf/kBT + k 2 ) a (describing conventional FPE and an 
FFPE in an infinite homogeneous system respectively) 
can be chosen to be the same. Exemplarily, exp(ikx) is 
the eigenfunction of free motion in both cases; we denote 
its eigenvalues by A^ p and A F FP respectively. Thus, if 
Cfp has a (nondegenerate) zero eigenvalue, whose eigen- 
function corresponds to a stationary state, the same holds 
for C a . The stationary states of the systems described 
by FPE and by FFPE therefore coincide. For closed sys- 
tems (no currents at infinity), the stationary state is that 
of thermodynamic equilibrium and is given by the Boltz- 
mann distribution. This is a general property of each 
subordination process, since a state stationary in t is also 
stationary in n. 

The solution of FFPEs under the given initial and 
boundary conditions can be obtained by means of an 
eigenfunction expansion, as is generally the case for nor- 
mal and subdiffusive motion p| Jl^ , p^ , [l^| . If <p m {x) are 
the eigenfunctions of the Fokker-Planck operator, then 
the solution of FFPE can be expressed as 



P(x,t) 
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where $ m (i) are the corresponding temporal decay 
forms. Here the difference between sub- and superdif- 
fusive FFPE gets to be evident: in the subdiffusive case 
$m(t) are solutions of a fractional ordinary differential 
equation 
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(A n being real and negative). Hence the Q m (t) are 
Mittag-Leffler functions ■ On tne other hand the 

superdiffusive FFPE (being of first order in time) leads 
to 



dt 



$ m (t) = -K, a) (-A m ) a $ m (i) 
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corresponding to a simple exponential relaxation, 
®m(t) — exp(— Ki a \ \X m \ a t). Thus, in the case of a 
discrete spectrum and of real, negative X m the Levy- 
flight FFPE retains the exponential nature of the re- 
laxation to equilibrium, a behavior typical for normal 
FPE, so that only the corresponding relaxation times 
change. For example, the relaxation behavior of a par- 
ticle in a harmonic potential, f(x) = —72;, follows im- 
mediately from a standard solution of the FPE fTs) : 
The eigenfunctions can be expressed through those of 
the Schroedinger equation, and the spectrum consists of 
a zero eigenvalue, Ao = 0, and of equidistant negative 
eigenvalues, A„ = —(^//kpT)?!. Since the spectrum of 
a Fokker-Planck operator with a harmonic potential is 
discrete, the relaxation is multiexponential. The equi- 
librium state of such a system (the eigenfunction corre- 
sponding to Ao = 0) shows a Boltzmann distribution. 
The longest relaxation time is given by the first eigen- 
value, Ai = - 7 /fc B T, so that t = (k B T/-f) a /K {a) . 

Another interesting example corresponds to the mo- 
tion in the absence of a field of a particle in an interval 
with absorbing boundaries at x = ±Z. The eigenfunctions 
of the Fokker-Planck operator are now the trigonometric 
functions, <j> m (x) — cos [(m + \/2)itx/l], and the corre- 
sponding eigenvalues are A m = — [(m + l/2)7r//] 2 . The 
eigenvalues of C a are A m = —Kr a -\ [(m + l/2)ir/l] 2a , so 
that the overall relaxation again follows a multiexponen- 
tial pattern. The survival probability for a particle ini- 
tially situated at the middle of the interval, x = 0, is 
equal to 
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At longer times this decay tends to a simple exponen- 
tial with the characteristic time r = KZ^ {l/2it) 2a . Note 
that the ^-dependence of this characteristic time dif- 
fers from that encountered in normal diffusion, where 
t = K~ 1 (l/2ir) 2 . In the case a — 1/2 a sim- 
ple analytical expression holds at all times: P(t) — 



arctan 



exp 



( TX K (l/2) 

V 2 1 



see Eq. 5.2.4.8 of Ref. []2(J. 

Using a representation of Levy flights in terms of a 
subordination of random processes and following the 
Kramers-Moyal procedure we have derived a fractional 
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Fokker-Planck equation for Levy flights, ft has been 
shown that when the regular Fokker-Planck operator has 
a discrete spectrum (as happens under appropriate po- 
tentials or boundary conditions) anomalous transport re- 
sults in an exponential relaxation towards an equilibrium 
distribution. These properties of the Levy-flight FFPE 
are compared with earlier findings for subdiffusive FFPE. 
The latter is characterized by a non-exponential Mittag- 
LefHer relaxation. The equilibrium solution corresponds 
in both cases to the Boltzmann distribution, emphasizing 
that there is no need to modify conventional thermody- 
namics in order to obtain strange kinetics. 
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